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Abstract 

Short-hard 7-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron 
star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in 
globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger 
is expected to be short, r ~ 0.1 Gyr, implying that the redshift distribution of merger events should follow that 
of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation 
and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays 
comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed 
DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed 
redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due 
to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are 
dominated by primordial or dynamically formed DNSs. 

Subject headings: gamma rays: bursts — stars: binaries — stars: neutron — gravitational waves 



1. INTRODUCTION 

Observations of 7-ray bursts (GRBs) indicate that they 
divide into two classes (Kouveliotou et al. 1993). 
GRBs of one class are of relatively long (> 2 — 
200 sec) duration and have softer spectra. Long-soft 
GRBs occur in star formi ng gal axies with high redshift 
z (van Paradijs et al. I 1997h . a nd their association 
i n several cases with type Ibc SNe tGalama et al. 1998: 

EStanek et al. 2003:" Hiorth et al. 2003 J Malesani et al. 2004: 
Campanaetal. 2006) suggests that they are the resul t of 
core collapse SN exp losions of massive stars fWoosle v 1993:i 
fPaczviiski 1998: MacFadven & Wooslev 1999.) . The second 
class of GRBs have short duration (< 2 sec) and harder spec- 
tra. Afterglows of short-hard GRBs (SHBs) have only re- 
cently been observed (e.g., Gehrels et al.; 2005; Bloom et 
alOoOS; Berger et al. 2005; Fox et al. 2005), and this has led 
to the first identifications of SHB host galaxies. In contrast 
to long GRBs, SHBs were found to occur in at least some 
cases in elliptical galaxies with very low star formation rates 
(SFRs), of order < Q.IMq yr~^. It is therefore unlikely that 
the progenitors of SHBs are also massive stars, since these 
have very short ('^ few Myr) life times. 

The gravitational wave (GW) driven merger of a double 
neutron star (DNS) may lead to a SHB (Goodman 1986; 
Pacz yiiski 1986; Eichler et al. 1989; Narayan, Paczynski & 
Piran 1992). If the time-lag t between DNS formation and 
the merger is large (comparable to a Hubble time), it forms 
a natural explanation why SHBs occur in galaxies where the 
SFR is very low. 
We consider two formation mechanisms for DNSs. If two 
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massive stars are born as a binary system, a DNS may form 
after the super nova (SN) explosion of both components ("pri- 
mordial" DNSs). Population synthesis models of massive bi- 
naries yield typically short merger times (r « 0.1 Gyr), im- 
plying that the SHB z-distribution function (DF) clos ely fol- 
lows the SFR history (e.g.. Bloom, Sigurdsson & Pols 199?; 
Belczynski et al. 200L 2006). The kick-velocity of the DNS 
is low (< 50km s^^; Dewi, Podsiadlowski & Pols 2005), and 
the spatial distribution of SHB mergers should follow the light 
distribution of their host galaxy (Portegies Zwart & Yungel- 
sor 1998; Bloom et al j 1999i) . 

Another possibility is that at the moment of star formation 
the neutron stars (NSs) are not in the same binary system, 
but one of the NSs is in a binary with a low mass main se- 
quence (MS) star. In globular clusters (GCs), such binaries 
are likely to have an exchange interaction with a single neu- 
tron star (Sigurdsson & Phinney 1995; Efremov 2000; Grind- 
lay, Portegies Zwart & McMilla ril 200a hereafter GPZM06), 
and thus form DNSs. A significant fraction 30%) of all 
NS mergers in the Universe may stem from such dynamically 
formed systems (GPZM06). 

For exchange interactions to occur, the stellar density must 
be very high. The delay time for dynamically formed DNSs 
is therefore mainly determined by the time until core-collapse 
(CC) of GCs, which is typically comparable to the Hubble 
time (©. We show that the predicted z-DF of SHBs is dif- 
ferent for primordial and dynamically formed DNSs, so that 
future z observations may determine which formation chan- 
nel (if any) is dominant. To date, the few SHBs with z- 
detection favor dynamical formation, although possible biases 
and small number statistics make a final conclusion at this 
point impossible (^. We discuss our results in ^ 

2. DISTRIBUTION OF MERGER TIMES 

The delay time r is a sum of the time tec until the dynamical 
formation of a DNS during core-collapse, and the time tow 
until the DNS merges. Here we determine the resulting delay 
function {dp/dT)dyn of dynamically formed DNSs. 

GPZM06 performed scattering experiments with 
scatters and sigmaS in the Starlab environment 
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(Portegies Zwart et al. to determine the cross-section 

for the formation of NS binaries when a NS interacts with 
a (NS, MS) binary. We use the orbital parameters of the 
resulting binaries to determine the resulting DF of GW 
merger times dpcw /dtow, which was found to be well fitted 
by dpGW / dtQw oc ^q^^. Typical merger times tow are very 
short compared to the Hubble time. 

The DNS formation rate per GC can be estimated as F « 
4 Gyr~^ neVi{Npi-/20), where ri = 10^ pc^^ng is the number 
density of NS stars, v — lOkms"^?;! is the velocity disper- 
sion, and A'pi is the number of progenitor binaries containing 
one NS (GPZM06). Since F cx n, DNSs form when the GC is 
very dense, i.e., during the CC phase. The delay between the 
formation of the GC and the SHB is the sum of the delay time 
icw and the time tec between formation of the GC and CC. 

For the DF of CC times we make the following assump- 
tions: the CC times of the GCs in our Galaxy are representa- 
tive for the whole Universe; the relation between tec and the 
half mass relaxation time is only a function of concentration, 
and it is well approximated by the relation given by Quinlan 
the formation rate of GCs is proportional to the to- 
tal SFR. We (conservatively) neglect repeated phases of CC 
("gravothermal oscillations"; Sugimoto & Bettwieser l 1983t 
Makinc 1996), which would lead to even lower redshifts. 

We use the half-mass relaxation times given by Harris ( 
□"996.) to find the cumulative DF Pcc(< tec) for GCs with CC 
times smaller than tec (Fig. The time tec for a GC between 
formation and CC is somewhat uncertain. Results in the lit- 
erature for single mass systems without binaries agree ap- 
proximatel y (e.g. , Ouinlan I 199 6: Joshi, Nave & Rasio 20011: 
Baumgardt i 20011) . A spectrum of masses can significantly 
decrease tec (Giirkan, Freitag, & Rasio 2004), while primor- 
dial binaries increase tec (Fregeau et al. 20031) . The estimate 
by Quinlan fl996) leads to most GCs having tec > tn (Fig. 

consistent with the observation that most (~ 80%) GCs in 
the Galaxy have not yet experienced CC. 

The probability function of the total time r = tew + tec to 
be in the interval (t, r + dr) is given by 



drJdyn 
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dt, 



dpc 
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atec Jo atcw 

The resulting delay function is shown in figure (|3; we 
weighed the CC times by the GC luminosity. Since the GW 
inspiral time DF diverges towards small times, {dp/dT)dyn is 
mostly determined by the DF of of CC times. It is also shown 
in fig. (|2j that the delay time DF {dp/dT)dyn is not modi- 
fied significantly by choosing tec which is smaller than that of 
fig-fflby a factor of 10. This demonstrates the robustness of 
the conclusion that for dynamical DNSs, the delay times are 
comparable to the Hubble time. 

In contrast to the delay function for primordial DNSs, we 
find that the delay function of dynamically formed DNSs 
grows with the delay time. For delay times shorter than 10 
Gyr, the average delay time is f 6 Gyr. 

3. REDSHIFT DISTRIBUTION 

We derive the predicted z-DF of SHBs, which depends on 
the event rate of SHBs as a function of z, the luminosity DF 
of SHBs, the delay function, and the detection threshold. 

3.1. The intrinsic z-distribution 

For DNS mergers, the intrinsic (as opposed to observed) 
SHB rate is given by the convolution of the star formation 
rate SFR(z) with the distribution dp(r)/dr of time delays. 
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Fig. 1 . — Cumulative DF P(< tec) of the CC times of the GCs in our 
Galaxy. Most GCs have CC times that are larger than 1 Gyr 
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Fig. 2. — The distribution {dp/dT)^iy^ of inspiral times (black solid line) 
from the dynamically formed DNS mergers, combining the delays from CC 
and GW inspiral. The average delay time is f Ri 6 Gyr. For comparison, we 
also plot the DF that would result in the (unlikely) case that tec would be 10 
times smaller than our estimate (red dashed line). Since the delay function is 
still a growing function of tec and has a large average value (t fa 5 Gyr), 
the conclusions on the resulting z-DF would not be affected, stressing the 
robustness of our results. 



Mntr(^) OC dz'-^SFRiz')^ [t{z) ~ i(z')] (2) 

(e.g. Guetta & Piran f 20051 [hereafter GP05] ; Nakar, Gal-Yam 
& Fox 200^; t{z) is the age of the Universe as function of z. 
We employ the SF2 model of Porciani & Madau (2001) 



SFRpm(2) oc 



(e3-4^ + 22)(l + z)3/2 



and the Rowan-Robinson \ 1999ft SFR, 

SFRrr(z) oc <^ ^ -^i^ 



(3) 



(4) 
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Table 1 

Model parameters and goodness of ht 



Fig. 3. — Intrinsic SHB rates as a function of z. For primordial binaries, 
tlie SHB 2-DF (dotted and dash-dotted lines) is only slightly shifted to lower 
z, compared to the SFR distribution (solid hnes), by the delay function of 
primordial DNSs. SHBs from mergers of dynamically formed DNS occur at 
much lower red-shift due to the large delay (dashed Hne). 



In figure we show the intrinsic SHB z-DF N(z) for 
(dp/(iT)pi.i„i oc 1/r (as possibly appropriate for primordial 
DNSs, Bloom et al. 1999; Belczynski et al. 120061) and 
dp/dr = (dp/dr)dyn (Eq. Q). 

3.2. The SHB luminosity function 

For the observed z-DF, the DF of peak luminosities and the 
minimal observable flux Pii,„ are required. Here we follow 



the procedure outlined in Guetta, Piran & Waxman £2004). 
We consider a broken power law peak luminosity function 
(LF) with lower and upper limits, 1/ Ai and A2, respectively: 



<^o{L)d\ogL = Corf log L 



(L/L*)-- 



L*/Ai<L<L* 
L*<L<A2L* 



(5) 

where Cq is a normalization constant. This is the "isotropic- 
equivalent" LF, i.e. it does not in clude a correction factor due 
to beaming. Following Schmidt f2001) we approximate the 
effective spectral index in the observed range of 20 or 50keV 
to 300 keV as -1.1 (N{E) cx E'^-^ and we use Ai,2 = 
(30, 100) (GP05). Both values are chosen such that even if 
there are bursts less luminous than L*/Ai or more luminous 
than A2L* they will be only very few (less than about 1%) 
of the observed bursts outside the range (L*/Ai, L*A2). We 
find a — 0.6 and P — 2, and then constrain L* by comparing 
the predicted peak flux distribution with the one observed by 
BATSE (see Guetta et al. [2004] and GP05 for details). The 
best fit values of i* are reported in Table Q. 

3.3. The observed z distribution 
The expected z-DF of the observed bursts is 

NinUz)dViz) f^— 



iVobs(^) 



1 + z dz 



<^o{L)d\ogL 



(6) 

where L^ax = A2L* — lOOi* and Lmin is the luminosity at 
z corresponding to the minimum peak flux Pn^^ required for 
detection. We estimate that for SWIFT is similar to that 
for BATSE, ^ 1 phcm^^ s^^, based on the observation 
that the detection rate of GRBs by SWIFT (- 100 yr"!) is 
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' The redshifts used are z={0.225, 0.16, 0.257, 0.55} for SHBs 
{050509, 050709, 050724, 051221). We consider separately the case 
including GRB050813, which is tentatively associated with a cluster 
at = 0.7 (Gladders et al. 2005). 



similar to that of BATSE (taking into account the different 
field of views and the fact that BATSE was triggering only 1/3 
of the time), and that the fraction of SHBs is similar for both 
BATSE and SWIFT. Fig. (0} shows a comparison between the 
observed and the expected integrated z-DF of SHBs for the 
different models listed in table (QJ. Our model for dynamical 
DNS mergers fits the observed SHBs much better than the 
model for primordial DNSs (see the results of Kolmogorov- 
Smirnoff [KS] tests in table (TJ). 

We note, however, that ruling out a primordial binary z- 
DF may be premature based on current data. For SWIFT, 
only 1/3 of detected SHBs have secure z determinations. If 
there is a bias against obtaining a secure redshift for higher z 
bursts, the observed distribution would be shifted, compared 
to the expected distribution, to low z. Assuming, for exam- 
ple, that obtaining a secure z requires a higher com- 
pared to that required for detection, the lower rate of de- 
tection of SHBs with secure z is accounted for by choos- 
ing Plim — 2.5phcm~^s~^ (which reduces the detection 
rate by a factor of 3 compared to that obtained for P\i^ = 
Iphcm^^s^^). Indeed, inspection of the 15 to 150 keV 
peak fluxes of SWIFT SHBs shows that the average peak 
flux of SHBs with z is higher than that of SHBs without z 
((/^) = 5.7 ± 5.7phcm-2s-i for GRBs 050509, 050724 
and 051221 compared to {fnoz) ~ 1-2 ± 0.5phcm^^s~^ 
for GRBs 050911, 051105, 051114, 05121 and 051227; 
Barthelmy et al. ■ 2005^ G ehrels et al. I 2005 ; Bloom et al. 
2005; Berger et al. l 2005t Fox et al. l2005t Villasenor et al. 
2005). 

Choosing Pii,„ = 2.5phcm^^s^^, the expected z-DF of 
primordial binary mergers is marginally compatible with ob- 
servations (see tableland fig. 1^. The time delay distribu- 
tion {dp/dT)dyn expected for dynamical DNSs yields a z-DF 
compatible with observations for both choices of Pum- 

4. SUMMARY AND DISCUSSION 

We have shown (fig |3} that the z-DF of SHBs expected 
from dynamical formation and subsequent merger of DNSs is 
markedly different from that expected from primordial DNS 
mergers (assuming {dp/dT)p,-i^ oc l/r). The large time for 
core collapse shifts the DF of dynamically formed DNS merg- 
ers to low z. The observed z-DF of SHBs strongly favors that 
expected for mergers of dynamically formed DNSs, as com- 
pared to that expected for primordial DNS mergers (fig 0] and 
table 1; see also GP05 and Nakar et al. 2006). 

However, current data do not allow to rule out a z-DF con- 
sistent with that expected for primordial DNS mergers, since 
redshifts were obtained only for a minority of the detected 
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Fig. 4. — Observed (histogram) and theoretical cumulative 2-DFs for dif- 
ferent SFR evolution (SF2 and RR) and different merger delay times (l/r 
expected for primordial DNSs, and (dp / dr) ^y^-^ expected for dynamical 
DNSs). Model distributions are shown for peak flux detection threshold 
Pum = Pi = lphcm~^s~^, which reproduces the observed detection 
rate of SHBs, and for Pnm = f 2 = 2.5phcm~'^ s~^, which reproduces 
the observed detection rate of SHBs with z determination . 



SHBs. This may be due to a bias against obtaining redshift 
information for high z (faint) SHBs (figE] table 1 and discus- 
sion at the end of §3.3). Future observations should allow to 
better constrain the z-DF of SHBs, and thus to differentiate 
between models. For example, detection of only a few high 



redshift (z > 2) SHBs would severely constrain the contribu- 
tion of dynamically formed DNSs. 

If the formation rate of primordial and dynamically formed 
DNSs are comparable (GPZM06), there will be an anti- 
correlation between z and off-set of the SHB from the center 
of its host galaxy, because DNSs formed in GCs will be closer 
in z, but farther away from their host center, since GCs reside 
in the halos of galaxies. The large time delay also implies that 
more DNS mergers will be observed by LIGO and VIRGO 
(Nakar et al. 2006) 

An alternative method for constraining the progenitors of 
SHBs is to consider their demography (Gal- Yam et al 200^ 
Zheng & Ramirez-Ruiz 2006). Zheng & Ramirez-Ruiz ( 
[2006) showed that in order to account for the preponderance 
of SHBs in elliptical galaxies, the delay times should be large, 
following dp/dr oc with S ~ 1.5. This method, which 
does not rely on the observed z-DF, also indicates delay times 
longer than expected for primordial binaries. We note that the 
specific frequency of GCs (number of GCs per unit luminosity 
in the V-band) is larger for elliptical galaxies than for spirals 
(e.g. Harri s. 1991.) by an order of magnitude. 
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